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Abstract 


The  bathythermo grams  in  the  subpolar  water  region  (inshore  side  of  the 
Gulf  Stream)  show  often  temperature  inversions,  while  those  in  the  Gulf  Stream 
and  the  Sargasso  Sea  do  not  show  such  features*  A  similar  situation  was  found 
in  the  Kuroshio  area*  The  T-S  curve  obtained  with  a  bathysalinograph  and  a 
bathythermograph  indicates  that  the  temperature  inversions  correspond  to  high 
salinity  and  thus  to  the  unstable  density  stratification.  Thus  the  inversions 
seem  to  be  caused  by  sinking  of  the  warm  and  saline  subtropical  (Gulf  Stream) 
water*  Such  water  is  driven  into  the  subpolar  water  region,  is  cooled  by  mixing 
or  evaporation  but  retains  its  high  salinity  and  sinks  to  an  equilibrium  depth# 

A  simple  mathematical  model  is  developed#  In  this  model  an  isolated  water  mass 
is  assumed  to  be  enclosed  in  a  small  balloon-like  membrane.  The  initial  density 
inside  the  balloon  is  different  from  the  one  outside*  It  is  assumed  that  tem¬ 
perature  and  salinity  of  the  inside  water  mix  with  those  of  the  surrounding 
water  but  mixing  of  temperature  is  faster  than  salinity.  This  model  can  explain 
not  only  sinking  of  the  subtropical  water  but  also  upwelling  of  the  isolated 
subpolar  water  in  the  subtropical  water  region  where  the  temperature  inversion 


is  not  found. 


1 


* 


1.  Introduction 


A  conspicuous  difference  in  the  water  structure  between  the  subtropical 
water  region,  including  the  Gulf  Stream  and  the  Sargasso  Sea,  and  the  sub¬ 
arctic  water  region,  is  seen  in  bathythermograms  obtained  in  these  two  regions. 

In  the  EASTWARD  cruise  of  July  ll-ll*,  1966  the  bathythermograms  at  two  neigh¬ 
boring  stations  (Stations  16  and  17)  with  8  nautical  miles  distance,  show 
such  difference  (Figure  l).  The  one  of  the  Gulf  Stream  (Station  16)  has  an 
upper  homogenous  layer  from  the  surface  to  about  20  meters  and  the  thermocline 
below  it,  while  the  one  north  of  the  Gulf  Stream  (Station  17)  shows  a  sharp 
thermocline  near  the  surface  and  two  temperature  inversions  below  it  (Figure  2A). 
In  order  to  study  the  vertical  distribution  of  salinity  corresponding  to  such 
temperature  inversions,  a  bathysalinograph  (on  loan  from  Navy  Mine  Defense  Lab, 
Panama  City,  Florida)  was  used  with  an  ordinary  BT  at  Station  19  which  is  also 
located  in  the  subarctic  water  region.  The  vertical  distributions  of  salinity 
and  temperature  at  this  station  are  shown  in  Figure  2B. 

Temperature- salinity  relationships  from  BS  and  BT  records  are  plotted 
in  Figure  3,  indicating  also  depths  in  feet  and  isopycnals.  The  T-S  curve 
shows  that  the  vertical  stability  is  positive  from  the  surface  to  i|0  ft.  and 
neutral  in  general  between  hO  and  120  ft.  but  has  a  large  negative  value 
between  55  and  58  ft.  mostly  due  to  decrease  of  salinity  downwards  rather  than 
the  temperature  inversion  at  this  layer.  Below  120  ft.  the  water  stratification 
is  generally  stable  but  again  the  decrease  of  salinity  causes  large  instability 
between  130  and  liiO  ft.,  though  the  temDerature  inversion  may  also  contribute 
to  such  instability.  Frequent  presence  of  temperature  and  density  inversions 
in  the  upper  layer  in  the  subarctic  water  (the  Oyashio  Current)  in  the  North 
Pacific  Ocean  was  reported  by  Japanese  Oceanographers  (Kawai,  1955;  Kuroda, 

1958).  They  also  noted  that  the  inversions  were  not  observed  in  the  sub- 
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tropical  water  region  (the  Kuroshio  Current).  A  glance  of  the  T-S  curve  of 
the  subpolar  water  such  as  Figure  3,  suggests  that  the  temperature  inversions 
is  caused  by  an  admixture  of  a  foreign  watermass  (the  subtropical  water). 

Such  admixture  may  occur  also  in  the  subtropical  water ,  as  postulated  by 
Rossby  (1936)  who  considered  that  the  subpolar  water  is  transferred  by  lateral 
mixing  through  the  Gulf  Stream  into  the  Sargasso  Sea.  However,  there  is  no 
definite  evidence  of  presence  of  the  temperature  inversions  in  the  Gulf  Stream 
and  the  Sargasso  Sea. 

Generation  of  temperature  inversions 

It  was  often  observed  that  eddies  with  horizontal  dimensions  of  tens  to 
hundreds  of  kilometers  are  present  near  the  shear  zone  along  the  edge  of  a 
strong  ocean  current  (Ichiye,  196f>).  The  eddies  in  the  northern  or  inshore 
side  of  the  western  boundary  current,  such  as  the  Kuroshio  and  the  Gulf  Stream, 
are  mostly  of  warm  water  and  anti-cyclonic,  while  those  in  the  southern  or 
offshore  side,  are  of  cold  water  and  cyclonic.  Particularly  warm  eddies  were 
frequently  located  in  the  area  such  as  off  Cape  Hatteras,  where  the  western 
current  separates  from  the  continental  shelf.  These  eddies  are  usually 
smaller  in  size  and  shallower  than  cold  cyclonic  eddies.  They  seem  to  have 
shorter  lifetime  since  the  horizontal  stability  is  smaller  due  to  their 
negative  ( an ticyc Ionic )  absolute  vorticity  (Ichiye,  1966).  On  the  other 
hand  cold  eddies  are  very  stable,  not  only  in  horizontal  but  also  in  vertical 
direction,  because  upward  displacement  of  the  isopycnals  near  the  center  of 
the  eddies  accompanying  with  the  cyclonic  circulation  increases  the  density 
gradient  below  the  upper  homogeneous  layer  (Ichiye,  1965). 

Common  occurrence  of  the  temperature  inversions  in  the  north  or  inshore 
side  of  the  western  boundary  current,  may  be  related  to  frequent  presence  and 
disappearance  of  warm  eddies  which  consist  of  the  subtropical  water.  The  sub- 
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tropical  water  of  the  upper  layer  of  the  Gulf  Stream  or  the  Kuroshio,  is  driven 
into  the  subpolar  water  area  as  a  warm  eddy  either  by  meander  of  these  currents 
or  by  wind  action*  The  subpolar  water  is  cooled  either  through  the  surface  or 
by  horizontal  mixing  with  the  new  environmental  water*  However,  salinity  is 
in  general  slower  in  mixing  and  diffusion,  both  of  molecular  and  turbulence 
scale,  than  temperature,  as  often  seen  in  a  fact  that  the  vertical  distributions 
of  salinity  preserve  finer  structure  than  those  of  temperature*  The  cooled 
subtropical  water  with  higher  salinity  than  the  surrounding  water  is  heavier 
and  must  sink*  Also  in  an  anticyclonic  eddy  the  geostrophic  effect  tends  to 
displace  isopycnals  downwards*  These  two  effects  cause  the  upper  water  to 
sink  to  the  depths  where  the  temperature  inversions  are  found* 

On  the  other  hand  if  the  subpolar  water  is  driven  into  the  Gulf  Stream, 
mixing  with  the  surrounding  water  may  increase  the  temperature  faster  than 
the  salinity  and  thus  decrease  the  density.  Also  if  entrainment  of  the  sub¬ 
polar  water  is  made  as  a  form  of  cyclonic  eddies,  the  geostrophic  effect 
causes  upward  displacements  of  isopycnals  below  the  upper  homogeneous  layer* 
Therefore,  in  the  subtropic  water  region  the  transported  subpolar  water 
usually  remains  in  the  upper  surface  and  is  diffused  by  mixing* 

3*  A  balloon  model  of  vertical  motion 

In  order  to  study  vertical  motion  of  a  mass  of  water  transported  into  an 
area  of  a  different  kind  of  water,  the  mass  is  considered  to  be  enclosed  in 
a  balloon*  It  is  also  assumed  that  its  membrane  is  so  flexible  and  permeable 
that  the  pressure  on  both  sides  is  the  same  and  heat  and  salinity  is  freely 
exchanged  with  the  surrounding  water  (Eckart,  I960)* 

The  balloon  having  been  at  the  level  z,  is  displaced  to  the  level  z  £  £* 

The  dens ity^^out side  the  balloon  changes  fromp(£)fcop^t  jr j When^P  is  small, 
the  variation  of  the  outside  density^  p is  given  by 

b?o  =g 


(i) 


v . 
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Density  of  the  sea  water  is  dependent  on  both  temperature  T  and  salinity  S 


and  thus  y  can  be  expressed  by 

'  c 

P0  ~  3^5  (  1  T°  +  f1  So  )  (2) 

where  y  is  a  standard  value  of  density  and  and  Bare  thermal  expansion 
coefficient  and  contraction  coefficient  for  salinity  respectively*  If  the 
balloon  moves  very  slowly,  change  of  the  inside  density  V-  is  partly  due  to 
the  change  in  pressure  and  partly  due  to  change  in  temperature  and  salinity 
by  mixing  with  the  outside  water*  Therefore,  we  have 


$yL  =  G  jb/c*  h  Jr5  (  S  Tl  r  p  S  S  l  )  (3) 

where  c  is  the  velocity  of  sound  in  the  water,  and  8  £>,  §  "Jr  and^S^are  change 
in  pressure,  temperature  and  salinity  of  the  inside  water,  respectively*  The 
pressure  change  is  given  by 


(h) 


The  change  in  temperature  or  salinity  is  given  by  a  solution  of  the 
transport  equation  with  proper  boundary  conditions  which  depend  on  the  shape 
of  the  balloon  and  on  the  distribution  of  temperature  or  salinity  of  the  outside 
water* 


Since  there  is  no  information  on  the  shape  of  the  balloon,  a  simplified 
equation  instead  of  the  exact  transport  equation  will  be  used  in  order  to 
represent  the  mixing  process*  One  such  equation  is 


St;  /H  1 — i n-  (Tl  -  r«) 

SS  l  /  St  =  "n  (Sl  -  So) 

where^/Jjyt  indicate  the  change  in  the  Lagrangian  method  of  description* 


(5a) 

(5b) 

The 


coefficients  m  and  n  are  in  general  different  even  when  the  conventional 
diffusion  coefficients  are  the  same  for  temperature  and  salinity,  because 
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distributions  of  these  properties  in  the  outside  water  are  different.  If  the 


vertical  velocity  w  of  the  balloon  is  assumed  to  be  constant  during  its  dis¬ 


placement  by  ,  the  outside  temperature  TQ  may  be  expressed  by 


z? 


To  Cz+g')  =  fo  (z)  t  S'Slo  /3Z-  =  i  o  ■+-  ^ 't  3Tq]Q\ 

Relation  (6)  may  be  substituted  in  of  equation  (£a).  The  solution  of 
(5a)  after  time^'t during  which  the  balloon  is  displaced  by^  is  given  by 


6) 


st=(t 


Bto  -TrH|-e~TnSt 
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(7) 


^  BZ.  ’  ~  ’  ’  ’  ""  '  '  ~"3 

where  T,-  ^  is  the  initial  value  of  T  •  •  The  solution  of  (5b)  is  similar  to 

(7)  with  n  in  place  of  m.  The  buoyant  force  at/thjis  proportional  to  the 

difference  of  9^  and  at  £  f ,  This  difference  is  expressed  by 

%  U  *  ? )  -  ?;  (  Z  +  £)  -  (Z)  ~  fL  ( T  J  f  S  5>0  -  S  ?/ 
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(8) 


where  /\  ,  Tr>  and  S0  indicate  the  difference,  outside  temperature  and 
salinity  at  the  initial  level  z,  respectively. 

Further  discussions  on  equation  (8)  depend  on  the  coefficients  m  and  n, 
the  values  of  which  are  not  known  by  theory  only*  Two  extreme  cases  are 
treated  here.  When  hi)yp<<  | and.  Ti^  f«  |j  equation  (8)  can  be  written 

?0(TZ)T-(€HO  =  A(y°'^0  ^  £(9?°/dz  W  C'Z) 

J  +  ^  "St  (To  -  Tel)  -p?5v\  s-t  ('S.-S  ;,)(9) 

The  second  parenthesized  term  of  the  r.h.s.  of  equation  (9)  is  an  ordinary 
buoyancy  term  due  to  adiabatic  process  (Eckart,  I960)  and  the  third  and  fourth 
terms  indicate  correction  due  to  the  effect  of  mixing  between  the  inside  and 


outside  of  the  balloon.  Particularly  if  m  :  n,  the  sum  of  the  third  and 
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fourth  terms  equals  —  /\  /  p. _ p\*  This  quantity  indicates 

clearly  the  compensation  to  the  initial  density  difference  &  |  O  _  p 


J 


a 


J  6 


\ 


by  the  effect  of  mixing* 

When  vnSt>)(and  j  *  equation  (8)  becomes 

y0(^z)  =$ ? — f 


c 


t^§f  -  c 


-a. 

(10) 


a  2.  '  mbt-  97,  /  u 

Equation  (10)  indicates  that  the  initial  temperature  and  salinity  of  the 
inside  water  do  not  affect  the  buoyancy  at  -f  P  •  Also  rather  large  values  of 
Y'yvVt  and  n  S~t  (larger  than  3  if  95$  accuracy  is  maintained  in  neglecting 

mot)  or  <2  n  St}  against  1)  are  in  the  denominator  and  thus  the  effect 
of  vertical  gradients  of  temperature  and  salinity  of  the  surrounding  water  is 
much  less  than  for  the  pure  adiabatic  process* 

If  the  inertial  reaction  of  the  balloon  is  neglected,  the  initial  difference 
of  density  between  outside  and  inside  the  balloon  may  be  balanced  by  its  vertical 
displacement  which  can  be  determined  by  equating  equation  (8)  to  zero.  When 


m  z  n,  equation  (8)  may  be  simplified  as 


So'zW-?i (m) = A (?0~?je 'n% Pit  II’  +<?? ?L  C -\n) 

which  is  valid  for  any  range  of  )>i S*t  •  This  equation  is  again  a  modified  form 
from  the  adiabatic  case  which  yields  the  first  and  second  parenthesized  terms 
of  the  r*h*s.  of  equation  (8)* 

The  effect  of  difference  in  mixing  coefficient  m  and  n  can  be  explained 
by  two  examples  which  may  occur  when  the  upper  layer  of  the  subtropical  or 
subpolar  water  enters  into  the  region  of  the  opposite  water  mass.  First,  it 
is  assumed  that  the  subtropical  (or  Gulf  Stream)  water  with  temperature  28°C 
and  salinity  35>  °/oo  enters  into  the  subpolar  water  region  having  temperature 
23°C  and  salinity  33  °/oo*  The  sigma  t  is  22*U2  for  the  former  and  22*Ui  for 
the  latter  and  thus  the  initial  density  difference is  almost  negligible , 
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When  the  z-axis  is  taken  positive  upwards,  —  0*36  (°C  per  m)  and 

3  /<?  7,  =  -0.0U  (S/oo  per  m)  in  the  upper  3>0  m  of  the  subpolar  water  region 

near  Cape  Hatter as  as  determined  from  the  July,  1966  cruise  of  the  EASTWARD 

.  ✓**« 

(Ichiye,  1966b)*  The  values  ofC^and^  in  the  range  of  temperature  between  23° 
and  28°C  and  of  salinity  between  33  °/oo  and  3£  °/oo  are  3*1  x  10"^  g  cm"*3 
per  °C  and  7*5  x  10“h  g  cm“3  per  opt  respectively*  The  effect  of  adiabatic 
process  may  be  neglected  in  the  following  calculation,  because  it  causes  only 
the  change  of  0*0h  in  sigma-t  values  for  vertical  displacement  of  10  m*  The 
displacement  at  which  the  balloon  filled  with  the  subtropical  water  reach 
equilibrium  can  be  determined  from  the  following  equation  (  in  meters)* 

0.3/  (  ~5  -  0,56  %  /  rn  Sir)  (  i  -  <5  ~  'Tn  0  "*0 

=  0.75  (-2  +  0.0 a-  §  /r\bt)  0-  e~^^)  (12) 

If  m  «  n,  the  value  of  ^  satisfying  equation  (12)  is  nearly  zero  as  expected* 
Since  salinity  seems  to  be  slower  in  mixing  than  temperature,  it  is  assumed 
that  =  2  and  nS't  -  0*5«  Then  equation  (n)  yields^  =  -10.6  (m)*  Thus 
the  original  balloon  sinks  about  10  meters*  On  the  other  hand,  in  the  Gulf  Stream 
(:j> To almost  vanishes  from  surface  to  30  m  and  equals  0*091  (°C  per  m)  from 
30  m  to  100  m  and  B3 0fcZ-  -O.OC57  (°/oo  per  m)  from  surface  to  30  m  and 
equals  -0*003U  from  30  m  to  100  m*  When  the  subpolar  water  with  the  same 
temperature  and  salinity  as  in  the  previous  example  is  driven  into  the 
Gulf  Stream  region,  the  equation  similar  to  (12)  yields  the  vertical  dis¬ 
placement  for  equilibrium  as  C  =  65  (m)  for  surface  to  30  m  depth  and  s 

16  m  for  30  ra  to  100  m*  Although  these  values  of  O  are  significant  only  in 
okeir  OfQ&ifc!  of  magnitude,  thy  calculation  shows  tiiG  upward  displacement  oi 
the  original  balloon*  These  examples  may  illustrate  one  of  the  reasons  of 
existence  or  non-existence  of  the  temperature  inversions  in  the  subpolar  or 
subtropical  water  region  respectively. 
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The  equation  of  motion  of  the  balloon  is  derived  from  balance  among  the 


inertial  reaction,  friction  and  buoyant  forces*  Thus  it  is  given  by 

^  "  §  9\  1  {  ?L  (Z±g)}  =  0(l3) 

where  the  term  in  the  brackets  is  given  by  relation  (8)  with  t  replacing  cyr  * 

If  m  s  n,  the  bracketed  term  becomes  simple  as  indicated  before  and  equation 

(13)  becomes 


d2§ 


r  r  \ 


4~  ~  3  I  p  :  « 

d  t 


cL  t* 

When  mt«0,  equation  (lb)  becomes 


1  ?c  4  ?cdZ 


<W 


-d-  4- 1/ 
d±*-  *  d-t 


where 


,  +N4  =iL  *(5>o-?0 


(15) 


N  *  «-9( 


~L  ~d  9c 

it 


a  z. 


t  ^  c-2) 


(16) 


The  quantity  N  is  the  Vaisala  frequency*  Equation  (lit)  expresses  motion  of 
the  balloon  caused  by  an  initial  impulse  given  by  its  r*h*s.  When  ■yyi't  V>  1, 
equation  (lit)  becomes 


d  t t 


4-2.  f  k 


il  -L  f-j?l  + 
d-t  t  {  s-t-  T 


a 

-t 


l-)5  =  o 


(17) 


where 


CL 


a? 


O 


^  ^  3,  /L 

Solutions  of  equation  (1?)  can  be  expressed  by 

Kt/z 


(18a) 


(18b) 


_  A  M  K  ;  ,/z  ( s  t)  t. 6  Lvk,  G*))  (19) 

where 

are  numerical  constants  determined  from  the  initial  consition*  The  constants  x 
Send  s  are  positive  and  defined  by 


a-  — 

^  l  /  4  '  ,K  '  VjL  .  .  ■  _  - 

M  and  W  -_are  Whittaker Ts  function  (Erdelyi,  195#)  and  A  and  B 

'  (x,  X  2 


s4  =  B  4-  Kz 


(20a) 


•  r  * 
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(20b) 


When  st  >>  1,  an  asymptotic  expansion  of  (p>t)  includes  a  term 

/  V  -  , 

Therefore,  the  solution  (19)  includes  a  term 


u. 


{  -  Y  -t  +  J  -tj  <21> 

which  becomes  infinite  as  increasing  t.  Therefore,  the  motion  represented  by 

(17)  is  ultimately  unstable.  This  can  be  seen  directly  from  equation  (17) 

£ 

***— ^ 

which  yields  -b^/b  for  'the  coefficients  of^as  t  increases  infinitely.  However, 

if  becomes  large  enough,  the  equation  (lb)  is  not  valid  because  this 
equation  is  derived  with  an  assumption  of  small  amplitudes  of  ^  • 

Concluding  remarks 

More  rigorous  analysis  on  generation  of  the  temperature  inversions  may  be 
done  by  Eulerian  method.  It  leads  almost  inevitably  to  a  nonlinear  problem. 

Further  development  for  studying  this  problem  may  be  hinged  upon  deter¬ 
mining  the  process  of  dissipation  or  decay  of  eddies,  particularly  warm  and  anti- 
cyclonic  geostrophic  eddies.  Several  approaches  may  be  suggested.  A  person 
with  a  great  amount  of  money  and  a  little  bit  of  valor  might  use  an  instrumented 
submersible  to  track  such  an  eddy.  One  with  a  great  amount  of  valor  and  a 
little  bit  of  money  may  use  a  raft  which  floats  with  the  eddy.  One  with  a 
great  amount  of  money  and  no  valor  may  use  instrumented  neutrally  buoyant  floats. 
One  with  a  little  bit  of  money  and  no  valor  may  solve  the  nonlinear  problem 
with  a  high-speed  computer.  One  without  money  and  valor  might  enjoy  studying 
the  mathematical  problem  related  to  such  processes. 

This  study  is  supported  by  a  research  contract  AT  (30-1)  2663  with  the 
Atomic  Energy  Commission.  The  data  of  Figure  2  and  Figure  3  were  obtained  by 
the  E- 27-66  cruise  of  the  EASTWARD  supported  by  the  Duke  University  Marine 
Laboratory  Cooperative  Oceanographic  Program  through  NSF  Grant-G-17669.  The 
bathysalinograph  is  on  loan  from  U.  S.  Navy  Mine  Defense  Lab,  Panama  City,  Fla. 


-10- 


References 


Eckart,  C#  (I960)  Hydrodynamics  of  Oceans  and  Atmospheres  Pergamon  Press  pp*  290 

Erdelyi,  A*  (Editor)  (1953)  Higher  Transcendental  Functions  Vol.  I  McGraw-Hill 
Book  Co*  pp*  300 

Ichiye,  T.  (1965)  Geostrophic  Eddies  Part  I  Technical  Report  No*  CU-21-65  to 
the  Atomic  Energy  Commission  (unpublished  manuscript)  pp*  26 

Ichiye,  T.  (1966a)  Hydrographic  structure  of  an  eddy  on  the  edge  of  the 

Gulf  Stream  off  Cape  Hatteras  (Abstract)  Trans*  Amer*  Geophys*  Union  U7  (l) 

pp*  109 

Ichiye,  T.  (1966b)  Gulf  Stream  Bulletin  No*  5>  Lamont  Geological  Observatory 
(unpublished  manuscript)  pp*  7 

Kawai,  H*  (1955)  Notes  on  the  inversion  of  the  upper  water  temperature  in  the 
waters  to  the  northeast  of  Japan  in  summer,  Bull*  Tohoku  Reg*  Fish,  Res# 

Lab.  6,  71-80 

Kuroda,  R#  (1958)  Notes  on  the  inversion  of  water  density  in  the  vertical  dis¬ 
tributions  of  sea  water  density*  Bull.  Tokoku  Reg.  Fish  Res.  Lab*  11, 

82-87  (in  Japanese) 

Rossby,  C#-G*  (1936)  Dynamics  of  steady  ocean  currents  in  the  light  of 

experimental  fluid  mechanics.  Papers  in  Phys*  Oceanography  Met.  6  (l), 

PP*  h3 


*»  t 


M 


L 


-11- 


Explanation  of  Figures 


Figure  1  -  Track  chart  of*  the  EASTWARD  cruise  on  July  11th  to  July  15>th. 

Each  station  indicates  BT  station.  Hydrographic  stations  are  indicated 
with  The  four  digits  at  each  BT  station  indicate  the  time  when  BT 

was  hauled  up  and  the  figures  within  parenthesis  indicate  depth  in  meters* 
The  operations  at  hydrographic  stations  are  distinguished  with  letter  codes* 
"N"  is  for  Nansen  cast,  "P"  for  plankton  tow,  "V"  for  Van  Dorn  sampler  cast, 
nG"  for  gravity  core,  nIn  for  Isacs-Kidd  midwater  trawl,  "U”  for  underwater 
photograph  and  nCM  for  current  measurement  with  a  Roberts  current  meter. 

Figure  2  -  (A)  Examples  of  BT  traces* 

(B)  BS  (bathysalinograph)  and  BT  traces  at  station  19* 

Figure  3  -  T-S  diagram  at  Station  19  determined  from  BT  and  BS  records  shown 
in  Figure  2  (B). 
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